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Introduction: Much of the interest in the polar re- 
gions centers on the fact that they likely contain the 
best record of Martian climate change on time scales 
from years to eons. This expectation is based upon the 
observed occurrence of weathering product deposits 
and volatile reservoirs that are coupled to the climate. 
Interpretation and understanding of these records in- 
quires understanding of the mechanisms that involve 
the exchange of dust, water, and carbon dioxide be- 
tween the surface and atmosphere, and the atmo spheric 
redistribution of these species. We will summarize our 
use of the GFDL Mars general circulation model 
(MGCM), to exploration aspects of the interaction be- 
tween the global climate and the polar regions. For ex- 
ample, our studies [1] have shown that while the north- 
ern polar cap is the dominant seasonal source for water, 
it can act as a net annual source or sink for water, de- 
pending upon the cap temperatures and the bulk hu- 
midity of the atmosphere. This behavior regulates the 
annual and global average humidity of the atmosphere, 
as the cap acts as a sink if the atmo sphere is too wet 
and a source if it is too dry. We will then focus our 
presentation on the ability of the MGCM to simulate 
the observed diurnal variations of surface temperature. 
We are particularly interested in assessing the influ- 
ence of dust aerosol and water ice clouds on simulated 
surface temperature and the comparison with observa- 
tions. Surface thermal inertia and albedo are critical 
boundary inputs for MGCM simulations. Thermal iner- 
tia is also of intrinsic interest as it may be related to 
properties of the surface such as particle size and sur- 
face character. 

Model: The GFDL Mars general circulation model 
simulates the circulation of the Martian amosphere 
from the surface to roughly 90 km [2], The MGCM in- 
cludes parameterizations for radiative transfer associ- 
ated with C0 2 gas and for aerosols. An arbitrary num- 
ber of aerosol populations can be transported by the 
simulated circulation. Dust may be injected at the sur- 
face using a prescribed rate and spatial distribution. 
We have recently added a dust source scheme that 
associates injection with resolved wind stresses and 
parameterized dust devil activity. This scheme allows 
the seasonal cycle of air temperatures and dust to 
match observations well at times when large-scale dust 
storms are not occurring. The model has also proven 
capable of simulating global dust storms with interan- 
nual variability in size and timing of occurrence. A po- 
tential source of memory for interannual variability is 
the spatial distribution of dust on the surface, as sug- 


gested by spacecraft and telescopic observations of 
interannual albedo variations. An ongoing line of re- 
search is considering the coupling of hjection and 
sedimentation to the surface budgets of dust to inves- 
tigate their role in interannual variability and assess net 
transport of dust onto the polar caps. 

The water cycle is represented by surface ice and 
regolith water reservoirs, atmospheric transport and ice 
cloud formation [1,3]. The optical properties of pre- 
dicted ice clouds can be passed to the radiative heating 
codes, allowing cloud radiative feedbacks and dust- 
water ice interactions to be examined. 

Surface Temperature: 

The daily and seasonal variation in surface tempera- 
ture is a central element in the description of the mar- 
tian climate. In the case of an optically thin atmosphere, 
surface temperature provides the bottom boundary 
condition that fundamentally influences the profile of 
overlying atmospheric temperature. The low thermal 
inertia of the Mars surface allows for the large seasonal 
variation in diurnal-mean surface temperature that re- 
flects the seasonal migration of the subsolar latitude 
and the annual variation in insolation due to the eccen- 
tric orbit. We have used MGS TES surface tempera- 
tures and thermal inertia estimates [4,5] to derive ther- 
mal inertia and albedo maps suitable for use in the 
MGCM. It is important to note that estimates of thermal 
inertia must account for atmospheric opacity due to 
dust and water and C0 2 ice clouds. These effects are 
significant in the polar regions and will influence the 
characterization of the polar surfaces. By using rela- 
tively coarse spatial resolution compared to [4], we can 
more readily trade off spatial resolution for temporal 
resolution and relate the evolution of observed morn- 
ing and afternoon temperatures (and thermal inertia 
estimates) to variations in atmospheric opacity. 

Figure 1 shows the seasonal evolution of zonally- 
averaged daytime (2pm) surface temperature (con- 
toured) from a reference simulation representing rela- 
tively clear sky conditions. There is a large seasonal 
variation in temperature that reflects the seasonal ni- 
gration of the subsolar latitude and the annual variation 
in insolation due to the eccentric orbit. The advance 
and retreat of the polar C0 2 ice caps approximately fol- 
lows the 150 K isotherm. It is apparent that very strong 
temperature gradients develop along the retreating 
edge of the polar caps as spring advances into summer 
in each hemisphere. These gradients likely give rise to 
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strong local thermal wind systems that evidently are 
associated with observed local dust storm activity 
along the cap boundaries [6], 

Figure 1 also shows the difference between zonally- 
averaged TES surface temperatures and those from a 
reference MGCM simulation. This figure clearly high- 
lights the seasonal changes of observed temperatures 
that may largely be attributed to variations in atmos- 
pheric opacity. Temperature differences are minimal 
during the relatively clear NH spring/summer season 
when the atmospheric opacity assumed in the simula- 
tion most closely approximates that of the actual Mars 
atmosphere. The effects of a regional scale dust storm 
at L s =225° in the first mapping year and a major, planet- 
encircling dust storm at L,=185° in the second year are 
evident. A dusty dmosphere leads to an increase in 
morning temperature and a decrease in afternoon tem- 
perature . 

There are systematic temperature differences in the 
vicinity of the polar caps. These are due, in part, to 
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errors in simulating the polar cap latitude. Significant 
temperature differences are also due to the presence of 
dust and polar hood clouds in the vicinity of the polar 
caps. The aphelion season tropical water ice cloud has 
a clear influence on apparent tropical nighttime tem- 
peratures. We will show how simulated temperatures 
depend on atmospheric opacity. In a related manner, we 
consider how atmospheric opacity affects the determi- 
nation of surface thermal inertia. 
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Figure 1 . The seasonal evolution of zonally -averaged afternoon surface temperature anomaly derived from TES 
spectra. Afternoon temperatures nominally correspond to 1400 LT. Predicted surface temperatures from a MGCM 
simulation employing a low (x=0. 1 ) atmospheric dust column have been subtracted from the observed surface tem- 
peratures to highlight changes in surface temperature due to atmospheric opacity. The simulated surface temperature 
is contoured. A dusty and/or cloudy atmosphere leads to a decrease in observed afternoon temperature (and an in- 
crease in observed morning temperature) relative to the reference simulation. The effects of a regional scale dust 
storm in the first year (L s =. 225°) and a major, planet-encircling dust storm in the second year (L,=185°) are evident. 




